program test_my_tracking_code
use analyse_lattice
real(dp) fix(nvd)

call init(3,1,0,0)

call make_lattice
fix=0.d0
call find_orbit_with_fpp(lattice,1,fix)
write(6,*) fix


call make_twis_table_with_A(lattice,fix(3),"twiss_with_A.txt")
call make_twis_table_with_Invariant(lattice,fix(3),"twiss_with_Invariant.txt")


end program test_my_tracking_code
